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Abstract. Lightning in the solar nebula is considered to be one of the probable sources for producing the chondrules that 
are found in meteorites. Gamma-ray bursts (GRBs) provide a large flux of y-rays that Compton scatter and create a charge 
separation in the gas because the electrons are displaced from the positive ions. The electric field easily exceeds the breakdown 
value of a 1 V irr 1 over distances of order 0.1 AU. The energy in a giant lightning discharge exceeds a terrestrial lightning flash 
by a factor of ~ 10 12 . The predicted post-burst emission of y-rays from accretion into the newly formed black hole or spin- 
down of the magnetar is sufficiently intense to cause a lightning storm in the nebula that lasts for days and is more probable than 
the GRB because the radiation is beamed into a larger solid angle. The giant outbursts from nearby soft gamma-ray repeater 
sources (SGRs) are also capable of causing giant lightning discharges. The total amount of chondrules produced is in reasonable 
agreement with the observations of meteorites. Furthermore in the case of GRBs most chondrules were produced in a few major 
melting events by nearby GRBs and lightning occurred at effectively the same time over the whole nebula, and provide accurate 
time markers to the formation of chondrules and evolution of the solar nebula. This model provides a reasonable explanation 
for the delay between the formation of calcium aluminium inclusions (CAIs) and chondrules. 
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1. Introduction 

Chondrules are typically millimetre-sized stony spherules that 
constitute the major component of most chondritic meteorites 
that originate in the region between Mars and Jupiter and which 
fall to Earth. They appear to have crystallised rapidly from free 
floating molten or partially molten drops. The large amount of 
chondrules present in chondritic meteorites and their possible 
crucial role in planet formation, highlights the importance of 
identifying the mechanism responsible for their formation. The 
heat source responsible for melting chondrules remains uncer- 
tain and pro pos ed processes have been recently reviewed by 
and Rubin ( 2000). One of the earliest proposals for 
the formation of ch ondrules is lightning caused by turbulence 
in the solar nebula ftVhiPDlell96d l Cameronl 19661) . Lightning 
is a widesp read phenom enon that is not restricted to terrestrial 
rainclouds llUmanll987l) . Spacecraft have detected lightning on 
Jupiter and other planets. This model is not favoured for chon- 
drule production for a number of reasons including the limita- 
tion to the build up of a large electrostatic potential by the gas 
conductivity, the br eakdown potential of the gas and the avail- 
able energy source ( De sch & Cuzzil200 0: Eisen hour & Busecld 
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Weidenschillind Il997l) . It has also been pointed out that, if 
lightning were the mechanism responsible for chondrule for- 
mation, it would have to op erate on a large sca le comparable in 
size with the whole nebula JPilippet al.ll992l) . 

Almost all the proposed heat sources are local to the solar 
nebula. One exception is that the precursor grains near the sur- 
face of the nebula were melted when they effic iently absorbed 
X-ray s and y-rays from a nearby GRB (iMcBreen & Hanlon] 
1999;). This mechanism can produce a large amount of chon- 
drules in the nebula (~ 30 Earth masses) but it has a low prob- 
ability of occurrence (< 0.1%). X-ray melting of material has 
recently been demonstrated in the laboratory for the first time 
using a powerful synchrotron dDuggan et al.ll2002ll2003l) . The 
model of the X-ray melting of material by a GRB did not in- 
clude lightning caused by Compton scattering of y-rays by the 
gas in the nebula or the y-rays from post-burst emission. The 
effect of this new process is to induce giant lightning discharges 
over the whole nebula. 



2. Absorption of y-rays in the nebula 

The absorption of y-rays in a gas with solar abundance is domi- 
nated by molecular hydrogen and it is sufficiently accurate here 
to consider only the absorption by H2. Compton scattering is 



2 



B. McBreen et al.: GRBs and giant lightning discharges 



20 40 60 



Fig. 1. The path length of y-rays (a) and range of electrons (b) 
as a function of energy in molecu lar hydrogen . The y -ray data 
is taken from the compilation of | Cullen et all J 1 989l) and the 
electron data from lBerger & Seltzerl ( 1964Th 



the dominant process in H2 up to 70 MeV above which pair 
production and trident production overtake Compton scatter- 
ing and only 27% of y-rays with energy of 170 MeV and 5% 
at 1 GeV undergo Compton scattering. For y-rays with energy 
above a few MeV, the kinematics of Compton scattering are 
such that most of the energy is taken by the electron, which is 
scattered in the forward direction. An incoming pulse of y-rays 
in H2 is gradually transformed into electrons that move further 
into the nebula leaving a cloud of positive charge in its wake. 
Here we do not consider the scattered photons that also tend to 
move in the forward direction. The attenuation length of y-rays 
and range of electrons in H2 is given in Fig. 1; the range of the 
latter in H2 is less than y-rays with the same energy. As the en- 
ergy increases, the range of the electrons become an increasing 
fraction of the y-ray Compton attenuation length, thus facilitat- 
ing a larger charge separation. This p rocess has simila rities to 
nuclear explosions in the atmosphere dLongmirell 978). 

A charge separation will also occur with pair production 
because some positrons annihilate in flight creating a moving 
excess o f negative charge that leaves b ehind a positive excess in 
the gas dAskarv anil 19621 Ijellev et alJll966l) . The cross-section 
for positron annihilation is about 1 Barn/T where T = E/moc 2 . 
About 10% of the positrons with energy E = 400 MeV will 
annihilate in flight. In order to make progress with charge sep- 
aration it is necessary to adopt a model of the nebula. There is a 
reasonably well developed theory of the solar nebula including 
the formation of terrestrial and gas giant planets. The discov- 
eries of extrasolar planets and disks orbiting young stars are 
leading to the characterisation of protoplanetary disks around 
young stars. In most models of the solar nebula, after the early 
evolution of the Sun, accretion slows down and the nebula be- 
comes essentially quiescent. We adopt the power law relation- 
ship for the surface density S of the solar nebula as a function 
of radial distance r from the Sun given by 



S = S [r/L4t/]" 



(1) 



where So has the norm alising value of 4250 g cm 2 at 1 .0 AU 
and exponent n is -1.0 JCam eron 1995). The choice of So and 
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Fig. 2. The surface density, S, as a function of height above the 
midplane, H, in astronomical units (AUs) for two models with 
normalising values, £0 and exponent n, with So = 4250 g cm -2 
and n = -1.0 (solid line) and So = 1700 g cirT 2 and n = -1.5 
(dashed line)and for three radial distances of (a) 1 AU, (b) 3 
AU and (c) 5 AU. The diamonds and squares give path lengths 
of 150 and 300 g cirT 2 measured from the top of the nebula. 



n varie s between the models of the solar nebula. lHavashi et alJ 
ll 19851) proposed So = 1700 and n = -1.5 in the minimum-mass 
model of the solar nebula and these values a re consistent w ith 
minimum-mass models of extrasolar nebula jKuchneJ2004l) . 

In conventional models it is usually as sumed that lightning 
occurs in the dus ty midplane of the nebula ( iHoranvi et all 19951 
lLove et al 1 ll995l) . However a GRB will preferentially interact 
with the outer region of the nebula that is in the direction of 
the GRB source. Therefore it is necessary to model the vertical 
structure of the nebula to obtain the surface density perpendic- 
ular to the midplane. 

In the thin disk approximation the midplane isothermal 
pressure is given by JCameronl 19951) 



P c =0.25SQc s 



(2) 



where c s is the speed of sound and O =(GM/r)5 is the 
Keplerian angular velocity and M is the mass of the Sun. The 
vertical structure of the disk is given by 



p = p exp(-z 2 /H 2 ) 



(3) 



where z is the vertical distance from the midplane. po and H are 
the density at the midplane and scale height of the nebula and 
are given by po = 2.44 P c (NkT) -1 and H = n c s (2Q)~' respec- 
tively where N is Avogadros number and k is the Boltzmann 
constant. 

The surface density profile perpendicular to the midplane is 
plotted in Fig. 2 for r — 1, 3 and 5 AU using the two models. The 
surface density is reasonably constant in the midplane region 
and drops by only a factor of 2.7 in the first scale height and by 
the larger factor of 20 in the second scale height. Approaching 
the nebula from above the midplane, a path length of 300 g 
cm~ 2 will reach the midplane at 5 AU, and 0. 1 AU above the 
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midplane at 3 AU using the minimum mass model (Fig. 2c and 
Fig. 2b). 

3. Prompt GRB emission and charge separation 

The rate of GRBs when averaged over the Hubble volume is 
~4x 1(T 7 yr 1 galaxy" 1 JZhang & Meszarosl2003l) . The local 
GRB rate is smaller due to the drastic decrease in star forma- 
tion at low redshifts and a value of 0.25 x 10~ 7 yr~' galaxy -1 
is widely used at z = 0. We adopt a value of 10~ 7 yr~' galaxy -1 
for the GRB rate and ~ 10 7 years for the lifetime of a pro- 
toplanetary disk. Hence there is a probability near 100% that 
any protoplanetary system will be irradiated by a GRB. The 
isotropic equivalent energy of GRBs is in the range 5 x 10 51 
to greater then 10 54 erg. A GRB with an output of 10 52 erg 
will deliver ~ 10 6 erg cm -2 to a protoplanetary system that 
lies at a distance equal to a galactic radius of 10 kpc. The high 
energy emission extends to well above 100 MeV with a spec- 
trum compa tible with an ex tension of the electron synchrotron 
component foingus 11200 lh . An anomalous radiation compo- 
nent in the energ y band between ~ 1-2 00 MeV was discovered 
in GRB 941017 jGonzalez et al.l2003l) . This component varies 
independently of the prompt GRB emission from 50 keV - 1 
MeV and contains three times more energy. 

To model the charge separation and electric field in the neb- 
ula, we adopt a GRB that gives 10 6 erg cm -2 with 10% of the 
energy between 20 MeV and 100 MeV, which is assumed to be 
in 65 MeV y-rays that yield Compton scattered electrons with 
energy of 50 MeV. The incident flux of y-rays is « 10 13 photons 
irT 2 and is attenuated exponentially in H2 using a mass absorp- 
tion co-efficient of // = 1.33 x 10~ 2 cm 2 g _1 (Fig. 3a (i)). The 
positive charge also declines exponentially in the same way be- 
cause each Compton event creates a positive ion. The electrons 
travel on average an additional 10 g crrT 2 into the nebula in- 
cluding only the ionization loss (Fig. 3a (ii)). The path length 
is further reduced by a retarding electric field of 1 V rrT 1 which 
is typical of the breakdown value (Fig. 3a (iii)). The charge ex- 
cess is the difference between the two distributions (Fig. 3b). 
The net positive charge is confined to a layer of thickness 10 
g cirT 2 whereas the negative charge is distributed over a much 
wider range. The net positive or negative charge is 11% and 
has a value of q « 10~ 7 C vnT 2 . Two percent of the net nega- 
tive charge is beyond 300 g cm 2 and penetrates deeply into the 
nebula. 

Electrical breakdown occurs when a normally insulating 
gas suddenly becomes conducting in a strong electric field. The 
gas pressure is particularly important because lower pressure 
reduces the voltage necessary for the discharge and increases 
the width of the discharge channel. Breakdown occurs when 
the electric field is strong enough that a free electron accu- 
mulates ~1 eV of energy betwee n successive collisions with 
gas molecules JPilipp et alJll992h . This condition is given by 
eE s (ncr) _1 ~ 1 eV where E s is the discharge electric field, n 
is the number density of gas molecules and cr ~ 10" 15 cm -2 
is the collision cross section. E s has a value of ~ 2.4 xlO 6 
V nT 1 f or air and ~ 2.0 xlO 6 V rrT 1 for H2 at atmospheric 
pressure jLove et alll995t> . The measured value of E s in thun- 
derstorms is about a factor of 10 lower and this difference is 
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Fig. 3. (a) The number of y-rays that Compton scatter and and 
the number of positive ions as a function of path length in H2 
for an incident flux of 10 13 photons irT 2 (i). The number of 
electrons as a function of path length taking into account the 
ionization loss of 10 g cirT 2 (ii) and also including a retarding 
electric field of 1 V rrT 1 (iii). (b) The charge excess as a func- 
tion of path length for the ionization loss (solid line) and also 
including the retarding electric field (broken line). The model 
used is for r — 3 AU and So = 4250 g crrT 2 with n = -1.0. 



often attributed to energetic runaway electrons from cosmic 
rays or radionuclides that prematurely trigger the discharge 
dGurevich et a l. I bOOlh . The value of E s scales with pressure 
and has a value of 20 V rrT 1 to 1 V rrT 1 at pressures between 
10~ 5 and 5 x 10~ 7 of atmospheric pressure, typical of disks 
in planetary forming systems. E, could have a higher value i n 
dust loaded regions near the midplane jGibbard et alJll997h . 
The Compton scattered electrons produce ions and electrons 
that increase the conductivity of the gas above that caused by 
cosmic rays and radionuclides (Desch & Cuzzi 2000, Love et 
al. 1995). The large scale and rapid formation of the charge 
separation prevents significant discharging by the gas conduc- 
tivity. 

The charge separation described in the nebula is analogous 
to a capacitor. The basic equation for a parallel plate capacitor 
yields the potential difference V = qd(e )~' between the plates 
separated by a distance d. V attains a value of 10 14 V for q = 
10~ 7 C irT 2 and a representative value of d = 0. 1 AU. The elec- 
tric field is ~ 10 4 V irT 1 and greatly exceeds E 5 by a factor of 
at least 10 3 even for a GRB at a distance 10 kpc and ~ 40 for 
a similar GRB at the distance to the Large Magellanic Cloud. 
As the GRB interacts with the nebula, the charge separation 
creates a s trong electric field that exceeds E, and triggers the 
discharge jDwverl2003l). The duration of the discharge is esti- 
mated at 100 s (iPilipp et al.ll 19921) . The average current in the 
channel is ~ 10 11 A if half of the excess charge from an area ~ 
d 2 flows down the channel in 100 s. The width of the lightning 
channel depends on the gas pressure and has an estimated value 
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of ~ 10 5 cm ass uming it is limited to a few thousand electron 
mean free paths dPilipp et alll992l) . 

The total energy in y-rays, over an area comparable in size 
to the charge separation of 0.1 AU, gives an upper limit of 10 29 
erg to the energy dissipated in the channel and this exceeds a 
large terrestrial lightning flash by ~ 10 12 . We cannot exclude 
the possibility of repeated strikes over the lightning channel, a 
situation that is somewhat analogous to the stepped leaders and 
return strokes in terrestrial lightning. Furthermore there is the 
possibility tha t the lightning may fragment into many channels 
(IUmanlll987l) . The visible and ultraviolet radiation from the 
discharge heats and melts the precursor grains to form chon- 
drules out to a distance of ~ 10 9 cm from the discharge channel. 
The total amount of chondrules produced in the nebula within 
a radius of 5 AU is 5 x 10 20 g assuming a GRB with 10 6 erg 
cm" 2 has 10 5 erg cm" 2 in 65 MeV y-rays, 2 x 10 10 erg g" 1 to 
heat and melt the precursor dustballs and a n efficiency of 10" 2 
to convert the y-ray energy to chondrules Jjones et afll2000l) . 
The amount of chondrules produced is too small to account for 
the total mass of ~ 3 x 10 24 g in the asteroid belt and the value 
of ~ 10 27 g when the asteroid belt was 300 times larger than at 
present. The mass of chondrules can be increased by 3 x 10 4 
to 1.5 x 10 25 g for a GRB at a distance of 300 pc (x 10 3 ) with 
a higher isotropic luminosity of 10 53 erg (x 10) and an anoma- 
lous MeV component as observed in GRB 941017 (x 3). This 

model is not sufficient because there is evidence from com- 

■ II I 

pound chondrules (Wasson 1993) and compositional gradients 
in chondrules JWasson & Rubinll2003l) that they were melted 
several times requiring a repeating process. 

4. GRB post-burst emission and other sources 

The afterglow continues for days and weeks after the GRB 
when the relativistic blast wave decelerates by sweeping up 
the external medium. Anomalies in the afterglow of some 
GRBs including the detection of Ge V y-rays that lasted for 
over 90 minutes after GRB 940217 JHurlev et all 1 1994 can 
be accounted for by post-burst emission. After the GRB there 
could be a more extended period in which energy is injected 
into the remnant. This activity may be caused by continued 
drainage of matter left over from the explosion into the newly 
forme d black hole or th e spin-down of a super-pulsar or mag- 
netar jRees et al.ll2000l) . The y-ray post-burst emission has 
not been measured an d we must rely on model computations 
jRamirez-Ruiz lE004T) . The relativistic outflow from a super- 
pulsar interacts with photons from a binary companion or from 
the explosion. The emission is mostly in y-rays and estimated 
at 10 48 erg s" 1 immediately after the GRB that declines to 10 47 
erg s" 1 after 10 hours and 10 45 erg s" 1 after 100 hours. The 
emission is beamed into a large solid angle that is not well con- 
strained by the model and we estimate it covers 10" 1 of the sky. 
The prompt GRB fireball emission is collimated into a much 
sma ller angle that is e stimated on average to be 2 x 10" 3 of the 
sky jFrail et all200lb . The nebula is 50 times more likely to be 
irradiated by post-burst emission than by the GRB. 

The model of GRB emission is not unique and in the case 
of a structured jet, where the energy density per unit solid angle 
falls away from the axis, the emission is beamed into a larger 



solid angle by a factor of about 5 JZhang & Meszarosl2003l) . In 
this case the post-burst emission is ten times more likely to im- 
pinge on the nebula than the GRB. The post-burst emission of 
up to 10 52 erg in y-rays is sufficiently intense to cause a major 
lightning storm in the disk that lasts for days or even weeks and 
up to 50 such events may occur over the lifetime of the disk. 
These fifty events, occurring at a distance of 10 kpc, will pro- 
duce 5 x 10 23 g of chondrules assuming an efficiency of 10" 2 
for conversion of y-ray energy to chondrules. The largest mass 
of chondrules is produced by the post-burst emission from the 
nearest GRB; a burst 100 pc distant will produce 10 26 g within 
a radius of 5 AU in a disk. Strong dependence on the distance 
between the GRB and the protoplanetary disk results in a wide 
range in the mass of chondrules produced. In the above esti- 
mates we have assumed that the charge separation is not re- 
moved by the conductivity provided by the ions and electrons 
produced by the Compton scattered electrons, cosmic -rays and 
radionuclides. This will only be correct if the post burst emis- 
sion is in short duration outbursts like the GRB emission. The 
post-burst emission is more likely to be highly variable in the 
case of continued and sporadic accretion into the newly formed 
Kerr black hole (McBreen et al. 2002). 

Soft gamma-ray repeaters (SGRs) are also a possible source 
of lightning in the nebula. The energy that drives the giant 
flares (> 10 44 erg) such as the 1979 March 5 event from 
SNR N49 may be caused by a sudden large scale arrange- 
ment of the magnetic field which releases magnetic energy 
jThompson & Duncan! 1 19961) . The extreme possibility is that 
the entire dipole moment is de stroyed in a single event releas- 
ing ~ 10 47 erg JEichler1l2002h . The rate of such events could 
be as high as the rate of magnetar production ~ 10" 3 yr" 1 and 
~ 10 4 such events could occur over the lifetime of the neb- 
ula. About 10 to 100 of these SGR explosions should be close 
enough to generate lightning in the disk and each one produce 
a small amount of chondrules, but the overall contribution of 
SGRs is much smaller than that of GRBs. 

Other nearby and less variable y-ray sources such as 
quasars, microquasars and powerful x-ray binaries (Mirabel & 
Rodriguez 1999) might be expected to generate charge separa- 
tion and lightning. They require a much longer time of ~ 10 5 s 
to generate an electric field that can exceed E s , but the gas con- 
ductivity may limit the electric field to below this value. There 
is sufficient energy in these sources to produce large amounts 
of chondrules provided a method is found to limit the gas con- 
ductivity. The same limitation has been noted in other models 
of charge separation and lightning in the protoplanetary disk 
(Desch & Cuzzi 2000, Gibbard et al. 1997, Love et al. 1995). 

5. Discussion 

GRBs and other sources are capable of producing significant 
quantities of chondrules in protoplanetary systems anywhere 
in the disk of a galaxy. Stars form in giant molecular clouds 
of size about 100 pc and hence there is a higher than average 
probability it will contain massive stars that produce GRBs and 
SGRs and irradiate protoplanetary systems in the same cloud to 
form large amounts of chondrules. 
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Compound chondrules consist of a primary that solidi- 
fied fi rst and one or more secondaries attac hed to the pri- 
mary dWassonl fl993l IWasson & Rubinl l2003h . Sibling com- 
pound chondrules have very similar textures and composi- 
tions and most, perhaps all, consist of chondrules melted in 
the same event. If this event is identified with a nearby GRB, 
the chondrules would have been produced all across the disk 
and should provide a simultaneous time marker between sib- 
ling chondrules in different meteorites. If the chondrules were 
produced by GRBs, the differences between the composition of 
meteorites are due to compositionally segmented regions of the 
nebula. The independent compound chondrules probably were 
produced in two separate events and it may be possible to de- 
termine the order and frequency of the GRB events from the 
meteorite record. 

The magneto-rotational instability (MRI) provides an un- 
derstanding of radial m ixing and turbulence in the disk 
dBalbus & Hawlevll998l) . Turbulence can concentrate dust par- 
ticles of a particular size to spatial densities well ab ove their 
background values dCuzzi et alJl20o4 IWoodlll997l) . A large 
amount of chondrules can then be produced by lightning and 
provide an explanation for the large quantity of compound 
chondrules in meteorites. The magnetic field in the nebula must 
be well coupled to matter for MRI to be effective and this con- 
dition was satisfied close to the Sun and beyond 10 AU where 
cosmic ray ionization may suffice to maintain a significant 
amount of ionization. The magnetic field during chondrule for- 
mation is not well constrained by the meteoritic evidence and 
seems to have a value b etween 1 and 10 gauss dSugiura et alJ 
ll979tljones et al!2000h . A magnetic field of this magnitude is 
significant because it would channel Compton scattered elec- 
trons into regions of enhanced magnetic field, such as MRIs, 
and possibly cause lightning between regions in the disk which 
is analogous to cloud to cloud lightning in the Earth's atmo- 
sphere. 

The 26 A1/ 27 A1 ratio of CAIs in meteorites has a value of 5 x 
10 5 whereas in chondrules the rati o has a much smaller value 
bones et alJbOOCi ICamerorJll995l) . The simplest explanation 
is that there was an interval of at least a million years between 
the formation of CAIs and chondrules. Chondrules should have 
formed preferentially late in the development of the disk after 
the Sun had stripped away most of the gas and allowed the y- 
rays to penetrate closer to the midplane. This effect is more 
important in models of the disk with more than the minimum 
mass as shown in Fig. 2. The penetration of the disk by y-rays 
is easier beyond the snowline at about 5 AU because the disk 
is less massive. Chondrules produced beyond about 5 AU may 
have formed earlier and hence could contain 26 Al. 



G RBs in the galaxy affect the Earths atmospher e dThorsetj 

1 19951 iKurt & Zaidellll99d IScalo & Whe"eierll2002l) and these 
events mav be responsible for mass extinc tions dDar et alJ 
ll998llMelott et alJ2004HThomas et al. l2004 . Fortunately the 
rate of GRBs is very low in our galaxy. However a GRB can 
reveal protoplanetary systems in other galaxies by transient in- 
frared emission fro m the melting chondrules and optical emis- 
sion from the gas I Duggan et alJ l2003t iMcBreen & Hanlon! 
Il999h . These events may be detectable for hundreds of years 



after the GRB when the expanding shell or cones of radiation 
interacts with protoplanetary systems. 
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